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Originally, it was noticed in the context of periodic potentials PjJ |2] that transformed 
Hamiltonians under special Darboux transformations are displaced with respect to the initial 
potential by half a period. Later on, this effect was studied in more detail with respect to first 
order Darboux transformations (also called 1-SUSY transformations) (31 Hj- It was shown that 
the range of possible displacements produced by a Darboux transformation is not restricted to 
half a period, but can take on values in a set of (in general) complex values called "Darboux 
displacements" (or SUSY displacements). The potentials allowing such displacement are called 
translationally invariant under the Darboux transformation or, simply, 1-SUSY invariant. It 
was proven that a real-valued even potential function Vo(x) = Vo(— x) is 1-SUSY invariant if 
and only if it is of the form V (x) = 2p(x + u>'), where p{x) = jp(x, g 2) #3) is the standard 
Weierstrass elliptic function with real and imaginary half-periods uj and u' and invariants g 2 , g% 
[Sj (including its degenerate forms such as 2x~ 2 , —2 cosh -2 x, and 2 sinh -2 x). This result really 
means that the family of 1-SUSY invariant potentials is rather sparse. Moreover, it was noticed 
in [3] that a simple displacement appears to be a "frustrated case" of the Darboux method. 
Nevertheless, this property has led to mathematically nontrivial results and to unexpected link 
between the theory of elliptic functions and supersymmetric quantum mechanics [31 E] . From 
a physical point of view a very remarkable property of 1-SUSY invariant potentials is that 
the corresponding Bloch solutions can be found analytically for any value of the energy. If 
a linear combination of Bloch solutions is used as the transformation function for a simple 
SUSY transformation, it produces an exactly solvable potential with locally perturbed periodic 
structure. We believe these potentials could find applications in describing contact effects in 
crystals, or modelling crystals with embedded inclusions. 

In this paper we continue the investigation of SUSY invariant potentials, but at the level 
of second order Darboux transformations. The aim of this work is to study what happens if 
a potential is assumed to be 2-SUSY invariant, i.e., a second Darboux transformation results 
only in a displacement Vq(x) — > V 2 (x) = Vq(x + d). From the very beginning it is clear 
that this condition is weaker than 1-SUSY invariance, since any 1-SUSY invariant potential 
is obviously also 2-SUSY invariant. Actually, as a preliminary result shows jl], the family of 
2-SUSY invariant potentials is richer than that of 1-SUSY invariant ones; at least, it includes 
a class of 2-soliton potentials which are not 1-SUSY invariant. Indeed, as we show below, it is 
even much richer, since a simple Darboux transformation over such a potential gives another 
2-SUSY invariant potential which cannot be reduced to a displaced copy of the initial one. 
Another remarkable property of a 2-SUSY invariant potential is that it allows for an analytic 
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representation of both linearly independent solutions of the corresponding Schrodinger equation 
at any value of the energy (not necessarily from the spectral set). Using this property we are 
also able to generate new exactly solvable potentials with locally distorted periodic structure, 
as illustrated below. Our results lead us to hypothesise that the general Lame-Ince potential 
is invariant under an nth order Darboux transformation. 

This paper is organized as follows. In the next Section we review several properties of 
Darboux transformations and 1-SUSY invariant potentials. In Section 3 we prove necessary 
and sufficient conditions for a potential Vq(x) to be 2-SUSY invariant and point out some 
simple implications of this property. Section 4 has purely illustrative character. Here we apply 
the results of the previous section to a well-studied family of 2-soliton potentials. Section 5 is 
devoted to an analysis of the n = 2 case of the Lame-Ince potential Vq(x) = n(n + l)p(x + u/). 
It should be mentioned that this remarkable potential has attracted much attention from 
mathematicians as well as physicists. Without going into further details we refer the interested 
reader to the excellent recent papers jH] where a vast literature is summarized. We believe that 
our general results applied to the n = 2 case establish new properties of this remarkable 
potential. In particular, we show that it is 2-SUSY invariant and give a simple analytic 
representation of its linearly independent Bloch solutions in terms of the Weierstrass functions. 
Further we give an explicit formula for 1-parameter family of two-gap (i.e. having only two 
finite forbidden and allowed bands) potentials. General properties of such potentials have 
been studied by algebraic-geometrical methods (see e.g. (HI El). The most striking of our 
developments is that we use only elementary means and the well-known properties of Weierstrass 
functions presented in 0. 

2. Preliminaries 

Darboux transformations (also known as SUSY-QM transformations) have become an 
important tool for dealing with spectral problems associated with the Schrodinger equation, 
especially in one spatial dimension. Though the basic procedure is quite simple, the method 
has attracted increasing attention from mathematicians and physicists for more than a century 



(see e.g. jEHH]). 



One begins with a Hamiltonian 




d = d/dx , x G M, 



(2.1) 
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(in appropriate units) whose eigenspace is two dimensional for any eigenvalue a G C: 

hou(a, x) = au(a, x) , h u(a,x) = au(a,x) . (2-2) 

As it is well-known, imposition of boundary, summability or Bloch type constraints selects 
"physical" solutions ip(E,x) and leads to physical interpretation of the eigenvalue E as a 
spectral parameter: 

hoi/>(E, x)=E ip(E, x), E e R . (2.3) 

If ui = u(a±,x) is a real and nodeless solution of equation (|2.2j) . then the 1-SUSY partner 
Hamiltonian 

hi = -d 2 + Vtix), V x {x) = V (x) - 2(ln u x )" (2.4) 

(the prime denotes the derivative with respect to x) is isospectral with Hq. (For non-periodic 
Vq(x) the spectrum of hi may change in one point, but we ignore this for the moment). 
Furthermore, an intertwining operator L obeying L ho = hi L has the form 

L = -«9+(log m)' . (2.5) 

Thus, the eigenfunctions of hi for eigenvalue a ^ ai are v (a, x) = Lu(a, x), hiv(a, x) = av(a, x). 
This procedure can be repeated by starting with hi to produce a second isospectral Hamiltonian 
/i2 = — d 2 + V2(x), etc. It can be shown that in this way, for fixed ui and «2 = u(ci2,x), the 
transformed potential is expressed in terms of the Wronskian W(x) = W(ui,U2) = Uiu' 2 — u[u2 
as 

V 2 {x) = V {x) - 2[log W{x)\" . (2.6) 

For a 7^ ai,a 2 the eigenfunctions of h 2 are v(a,x) = W(ui, u 2 , u) W / ~ 1 (ui, m 2 ), h 2 v(a,x) = 
av(a,x) where u = u(a,x). It is easy to see that W(u x , u 2 , m 1j2 ) oc u 2) i so that for a = ai,a 2 , 
up to an inessential constant factor, this formula gives 

vi,2 = vi <2 (ai <2 ,x) = u 2) iW~ 1 (ui, u 2 ) , h 2 vi j2 = a 1)2 v lj2 ■ (2.7) 

We note that this is a particular case of a much more general result for a chain of N Darboux 
transformations |9j. Solutions of the initial Schrodinger equation Uj = itj(cij,x) are called 
transformation functions while their eigenvalues aj are known as factorization constants. We 
refer the reader to [9J for more details. 

A particularly interesting case is where Vq is periodic: 

V (x + T) = V (x). (2.8) 
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We recall that inside the eigenfunction space of each eigenvalue for such a Hamiltonian we can 
build up Bloch functions, i.e., two linearly independent solutions m ± of /i M± = such that 

w ± (a,x + T) = [/3(a)] ±1 w ± (a,x) , /3(a) G C . (2.9) 

When for a = E one has \f3(E)\ = 1, the -E-values belong to a spectral band and the 
corresponding ("physical") solutions are bounded, while for \/3(a)\ ^ 1 they are unbounded 
( "non-physical" ) and the values of a lay in forbidden bands (see e.g. 10J). By using appropriate 
Bloch functions as transformation functions in the Darboux algorithm one can construct 
periodic partner Hamiltonians PO [I], [H]. A linear combination of Bloch functions leads 
to a perturbed periodic structure which is asymptotically periodic. 

Let, for instance, the transformation functions be Bloch functions Ui(x) and u^ix) with 
corresponding factors (3\^ = P(ai,2) as defined in ()2.9|) . Then, according to ()2.7|) . Bloch solutions 
of the transformed equation corresponding to eigenvalues ai and a 2 are respectively 



where 



v^(x + T)=fc 1 Vi(x), vt(x + T)=(3 2 vl(x). (2.11) 

If under a Darboux transformation we have 

V l {x) = V Q {x + d) (2.12) 

then d is called a Darboux displacement and Vo is said to be 1-SUSY invariant. The first cases 
of Darboux displacements were found for periodic potentials PI |2] . Later it was proven [3] that 
an even function Vq{x) allows for a first order Darboux displacement if and only if it satisfies 
the nonlinear differential-difference equation 

2 

= const , (2-13) 



V (x) + V (x + d) - - 



Vq(x) + V£(x + d) 



2 [V (x) - V (x + d) 

which, up to a constant, is equivalent to the addition formula for the Weierstrass p function. 
This result means that the family of even 1-SUSY invariant potentials is restricted to the p 
function and its degenerate forms, such as one-soliton potentials. In the next Section we prove 
necessary and sufficient conditions for a potential Vq to admit a displacement under a second 
order Darboux transformation. The potentials possessing this property will be called 2-SUSY 
invariant. 



Second Order Darboux Displacements 

3. Second order Darboux displacements 



6 



Let, as before, uf(x) and u 2 (x) be two Bloch eigenfunctions of h which are chosen as 
transformation functions for a 2-SUSY transformation. Although we concentrate on the case 
of periodic potentials Vq(x + T) = Vq(x), our results have a more general character, as will 
be mentioned below. Now we assume that the potential V 2 , obtained according to (|2.fi|) . is 
displaced from the original one, i.e., 

V 2 (x) = V (x + d). (3.1) 

Thus, the differential equations corresponding to ho and h 2 have exactly the same solutions, 
but shifted by the displacement d. Hence, taking into account the property (j2.11|) we see that 

v^(x — d) oc Ui(x) , v 2 (x — d) oc u 2 (x) . (3.2) 

^From this and (J2.1(Jj) it follows readily that 

W( ) M 2~( x ) u t( x ) ^2 ^ 

u^(x + d) u 2 (x + d)' 

where Ci }2 are constants. It is convenient to use the notation 

x = x + d, &(x) = l/W(x), p(x) = -2[ln (3.4) 
Therefore, p'(x) = V (x) — Vq(x) and from (|3.3j) it follows that 

Ui(x) = c u 2 (x) , (3.5) 

where c is an inessential non-zero constant. After taking the second derivative of (|3.5|) and 
using (|2.2j) . one obtains 

[V G (x)-a x ]<S>(x)u 2 (x) = [V {x)-a 2 ]<t>(x)u 2 {x)+u 2 (x)<S>"(x)+2[u 2 (x) ] V(z) (3.6) 

Now, from (J3.4|) one has 

/ 2 

$' = --$, $" = -^$ + ^$, (3.7) 
2 2 4 K J 

and from ()3.6|) this yields 

pnu 2 (a;)] = ^ + 4 + — ^ • ( 3 - 8 ) 



Similarly 
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2p 4 p 

n -/~m/ p' p a x -a 2 , . 

pn«j (re) = — - - + . 3.11 

2p 4 p 

Remark 1 TTie formulas $3.t^i - (\3. iij) illustrate one of the most remarkable properties of a 
2-SUSY invariant potential. The four Bloch functions corresponding to the eigenvalues a% 
and <22 producing a displacement are expressible only in terms of the potential difference 
p'(x) = Vq(x) — Vq(x + d) (a known function) and its primitive (also known from Eq. (\3.19i) 
below). 

Later on for a particular case of the Weierstrass potential we will integrate these equations to 
get analytic expressions for the solutions of the Schrodinger equation. 

We turn next to finding necessary and sufficient conditions for Vq(x) to admit a displacement 
under a second order Darboux transformation (shortly second order Darboux displacement). 



Proposition 1 The Hamiltonian ho allows for a second order Darboux displacement if and 
only if the potential Vq(x) satisfies the following differential-difference equation: 

2 



T 



4 



T_ 

M 



constant . 



where 



and 



T = (S'/Df + 4S- 2D"/D, M = D — (S'/D)' . 



S = V (X) + V (X 



D = V {x) - V (x) . 



(3.12) 



(3.13) 



(3.14) 



Note that this formula resembles (2.13) for the first order displacement, but the functions T 
and M involved depend in a much more complex way on Vq(x). 

Proof. Since u" = [(u'/u)' + (u' /u) 2 ]u, the Schrodinger equation for is 

,1 2 



"2 



+ 



u 2 



+ a 2 - V (x) = 



which together with (J3.8)) results in 



2pp" - {p'f + -p 4 + 4(ai - a 2 ) 2 + 2p 2 (a x + a 2 - S) = 



(3.15) 



(3.16) 
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where S is defined in ([3.14)1 . Now we take the derivative of ([3.16)1 and, since p' = V (x)—V (x+d) 
is a known function, we can treat the result as a quadratic equation for p, whose solution is 




P=-±\l( ¥ ) -2 P y + 4S-4(a 1 + a 2 ). (3.17) 

Then taking the derivative of ([3.17)1 . and rearranging, we obtain the nonlinear differential- 
difference equation for Vq(x) 

2 

(3.18) 



4(ai + a 2 ) = T - ~ 



T" 



where T is defined by (13.13)1 . Keeping in mind that the eigenvalues a>\ and a 2 are independent 
of x and taking into consideration the definition of M (second formula in ([3.13)1 ). we see that 
the necessary condition is fulfilled. We shall show now that this condition is also sufficient. 

First, we notice that by using ([3.18)1 to eliminate a\ + a 2 in ([3.17)1 we find 

p = ^ +1 2 D-ls'/Dy (3 - 19) 

Here we have to choose only plus sign in ([3.17)1 since it is easy to see that the minus contradicts 
the condition p' = Vq(x) — Vq(x + d). If the right hand side of ([3.18)1 is constant, then from 
([3.19)1 we obtain p which, by retracing our steps, satisfies ([3.16)1 which is equivalent to the 
Schrodinger equation ([3.15)1 for u 2 . Therefore, from ([3.8)1 and ([3.9)1 . we have the logarithmic 
derivatives of u 2 and uf. These being known, we are able to calculate the second logarithmic 
derivative of their Wronskian 

[in WY = ( 0l - ^flf-i - ( "''r,, ) 2 = \D, (3.20) 

[H U l/ U 2)\ \[H U l/ U 2)\J 2 

which by ([2.6)1 means that ho possesses a second order Darboux displacement. □ 

Corollary 1 Any 2-SUSY invariant potential allowing a set of displacements d generates a 
one-parameter family of 2-SUSY invariant potentials (some of them may have singular points) 
obtained by a 1 -SUSY transformation with one of the four function uf, uf as transformation 
function. 

Proof. Let Vq(x) be 2-SUSY invariant. This means that there exists a second order Darboux 
transformation with transformation functions u± and u\ which results only in a displacement of 
the transformed potential V 2 (x) = Vo(x + d). Any second order transformation may be factored 
into two successive first order transformations. One can realize the first transformation with the 
transformation function to get an intermediate potential V\(x) which cannot be obtained by a 
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displacement of Vq(x) if the later is not 1-SUSY invariant and hence V\{x) is essentially different 
of Vq(x). We shall show that Vi(x) allows for the same second order Darboux displacement. In 
order to show this property take the logarithms of the first equality in ()3.3|) 

2\ogW(x) + 2\ogUi(x + d) = 21ogci + 2\ogu 2 (x) . 

After taking two derivatives and using ()2.6|) and (|3.1|) we obtain 

V (x) - 2[log«2 (x)]" = V (x + d)- 2[\oguY(x + d)]" . (3.21) 

The left hand side of this equality means that we realize a 1-SUSY transformation of the 
potential Vq with the transformation function u 2 {x) and from the right hand side we learn 
that a similar transformation is realized of the shifted version of the same potential with the 
transformation function u±(x + d). Since the Darboux transformations are invertible, we see 
that these transformed potentials may always be related to each other with the help of a 2-SUSY 
transformation for which Vq is an intermediate potential. Hence, the equality (|3.21|) means that 
after such a transformation we get only a displacement by the value d. If the function u 2 (x) 
has nodes, the 1-SUSY partner of Vq{x) will have singularities (poles). To displace V (x) by 
another value we have to choose a second pair of transformation functions and u 2 that give 
us another member of the family of 2-SUSY invariant potentials. □ 

Remark 2 No member of the family obtained above is a displaced copy of the initial potential 
if it is not 1-SUSY invariant. Moreover, few of the potentials are displaced copies of other 
members of the same family, which means that the whole set of 2-SUSY invariant potentials is 
quite rich. 

Indeed, if the initial potential is not 1-SUSY invariant it is not possible to displace it by a 
1-SUSY transformation. Next, any two different members of this family are related by a 2- 
SUSY transformation, where Vq is an intermediate potential. The factorization parameters of 
such a transformation are independent of each other, which means that once a representative 
of the family is fixed, V, the whole family is characterized by two independent (factorization) 
parameters. The potential V, being 2-SUSY invariant, also admits 2-SUSY transformations 
resulting only in a displacement, but these transformations are characterized by one parameter 
only. This follows from equations (J3.16j) and (J3.18j) which may be considered as defining one 
of the factorization constants (say a 2 ) as an implicit function of the other (ai). This property 
will also be illustrated in the next Section and used in the Section 5 to produce new exactly 
solvable periodic potentials. 



Second Order Darboux Displacements 10 

Remark 3 The results of this Section are valid not only for a periodic potential but for any 
potential satisfying the conditions of Proposition 1. 



This follows from the fact that the only place where we used the Bloch property of solutions 
of the Schrodinger equation (that is the periodicity of the potential Vq(x)) is in formula (J3.2I ). 
Thus, and uf can be not only Bloch solutions, but two pairs of linearly independent solutions 
of the Schrodinger equation with a non-periodic potential. 

A simple sufficient condition for (ETH1 is that T = (S'/D) 2 + AS - 2D" '/D be constant. For 
example, consider the 1-SUSY invariant potential Vq(x) = 2p(x — u'). From the differential 
equation for p, we obtain Vq = 3Vq — g 2 and therefore D" = 3Vq(x) — 3Vq'(x + d) = 3DS. 
From ([2.13)1 . the 1-SUSY invariance requires (S'/D) 2 = 2S + const. Therefore, in this case, 
T = (S'/D) 2 — 2S = const and, consequently, V is also 2-SUSY invariant, which of course was 
evident from the outset. 

Note that for even potentials, Vo(— x) = Vo(x), the expression in square brackets in ()3.12|) is 
not defined for x = —d/2. Using Mathematica ^2] we have found that in this case 

lim i r r i 2 _ lgrong + to 2 c - vm H ? + Y^th + 

x ^ d/2 4 [D - (S'/D)'\ K)'[6(y ') 3 + V»V£» - U V (iv) ] 

where all the functions in the right hand side are evaluated at x — d/2. 

4. Two-soliton potentials 

In this section we illustrate our previous developments for a well-studied two-soliton potential. 
In its most general form this potential is determined by four parameters: 

2(af — a^)(al sech 2 (a 1 x + (3i) + a\ csch 2 (a 2 x + /3 2 ) (A I) 

[aitanh(aix + (3i) — a 2 coth.(a 2 x + (3 2 )] 2 

The parameters a 2 > ot\ > define the positions of the discrete levels Eq^ = — while 
fix and (3 2 characterize isospectral deformations. It has been shown by direct calculation [3] 
that the symmetrical case, fix — fl 2 — 0, is 2-SUSY invariant. We note that exactly the same 
calculations can be carried out for the general case, indicating that for any fixed values of the 
parameters ax, 2 and fix,2 the potential ()4.1|1 is 2-SUSY invariant. This illustrates Remark 3. 
With evident modifications of the results of [3] we give the solutions uf and producing the 
second order Darboux displacement d of Vq: 

u l,2 = W(uxo,U20,U 3 o i 4o)W~ 1 (uxo,U2o), 
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where u w = cosh(aiX + Pi), u 20 = sinh(a 2 x + /3 2 ) are solutions of the free particle Schrodinger 
equation generating the potential (|4.1j) from the zero potential, and ^30,40 = exp(a 3!4 x), 
«3,4 > 0. When parameters 0:3 and vary in a certain domain, independently of each other, 
these functions produce an isospectral deformation of Vo, resulting in the same expression with 
the different values of f3\ and p 2 . To select from these transformations only ones leading to 
the shift Vq(x) — > V 2 (x) = V (x + d) one has to impose an additional restriction on a 3 and 
«4. The analysis in reference 0] shows that to get real displacements one has to change the 
factorization constants a.3 and so that they fall between the existing discrete levels of Vq\ 

Ct\ < OJ3 < «4 < C( 2 . 

Let us now characterize the one-parameter family mentioned in Corollary 1. For this purpose 
we have to realize a 1-SUSY transformation of the potential Vq by using one of uf 2 as the 
transformation function. For simplicity let us fix Pi = p 2 = and choose the function u~l to 
produce the potential V\. According to properties of Darboux transformations (see e.g. 0) 
this is equivalent to a third order transformation of the zero potential with transformation 
functions W10, u 2 q and M30 or simply to a chain of transformations where, for instance, the first 
transformation is realized with the function M30. For such a transformation the zero potential is 
not affected, but the functions Wio and u 2 q receive changes and, up to a constant factor, become 
u\q{x) — > v w (x) = cosh(«ix + 71), u 2 o(x) — > v 20 (x) = sinh(a 2 a; + 72), where 



After the second order transformation with the functions v 10 and v 2 q we get from the zero 
potential the same two-soliton potential (j4.1j) where /3's are replaced by 7's. This is just the 
1-parameter family of Corollary 1 where is the parameter. From the definition (|4.2jl we see 
that the parameters 71 and 72 are not independent of each other: tanli7i/ tanli7 2 = a\/a 2 . 
They vary so that the potential V± cannot be reduced to a displaced copy of the potential 
()4.1j) at Pi = j3 2 = 0, since the later condition requires a\/a 2 = 71/72, which leads to 
tanli7i/ tanh72 = 71/72- The later is possible only for 71 = 72 which contradicts to a>i 7^ a 2 . 
Just from this, the potential Vi, being of course a two-soliton potential, is essentially different 
from Vq, as pointed out in Remark 2. 

We will illustrate other properties derived in the previous section for the simplest case of 
this potential. We choose a,\ = 1 and a 2 = 2 to produce the well-known Poschl- Teller potential 
well 



71 = arctanh (ai/a^) , 72 = arctanh {a 2 /a^) . 



(4.2) 



Vq(x) = —6 sech 2 x 



(4.3) 



for which E< 



— 



-4 and E 



1 — 



-1. 
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From (J3.13|) we find the value of T, 

T = 4 [-l + csch 2 d-4sech 2 x-4sech 2 (d + x)] , (4.4) 

which after being substituted into Eq. ()3.18|) gives the x-independent quantity 

a\ + a 2 = — 5 — 3 csch 2 <i . (4.5) 

This means that the necessary and sufficient conditions for Vo (|4.3|) to be 2-SUSY invariant are 

fulfilled. Now from (J3.19j) we obtain p(x) = —6 [tanh x + cschd cosh x sech (x + d)]. From Eq. 
(13.16)) we deduce the difference of the factorization constants 



a i — a 2 = ±(3 cothd v 1 — 3csch 2 d ). (4.6) 

From here and (J4.5j) one gets one factorization constant as a function of the other: a% = 
\ (—6 — a 2 ± a/— 12a 2 — 3a 2 ). It is clear that this equation and (j4.5|) define the displacement d 
as a function of a 2 : 



-4 - a 2 + v3A/-a 2 (a 2 + 4) 
d = arccsch— — . (4.7) 

It follows from here that real displacements are possible only when both factorization constants 
lie between discrete levels, —4 < ai j2 < —1, and in the above formula for a\ the lower sign is 
chosen. This agrees completely with the previous result We see also that when a± — > —3, 
then a 2 — > —3, and we get the minimal value of the displacement c? m i n = arccsh (l/v^3). The 
other limit is ai — ► —4, which gives a 2 — > —1, so that the displacement increases indefinitely, 
d —>■ oo. In this respect we remark that for 1-soliton potentials the (first order Darboux) 
displacements range is always (0, oo). 

We are able now to find the right hand side of equation ()3.8|) . It is expressed in terms of 
elementary functions only and its primitive gives us a solution of the Schrodinger equation with 
the potential V : 



x sech x 



\/2(coth d + 2 tanh x) sinh d + V cosh 2d — 7 



u -(x) = cos h(x + rf) e -2^ cotha! +Vi-^sch d } [ 3 _ cosh2 d + cosh2a; + cosh(2x + 2d)} 1 

1 1/2 

(4.* 



-y / 2(coth d + 2 tanh x) sinh d + V cosh 2d — 7 



The derivative of this expression yields the potential difference as a function of the parameter 
d: 



AV = [(4 cosh x csch cisech (x + d) + tanh x) (sech 2 (x + d) tanh(x + d) + 2sech 2 xtanhx) 
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+ (sech 2 (x + d) — sech 2 x) (3csch 2 d + coth d \J\ — 3csch 2 d (4.9) 

— sech 2 x — 2sech 2 (x + d) +3 tanh 2 x)\ [tanhx + csch d cosh x sech (x + cf)]~ 2 • 

Of course, this intermediate quantity can be eliminated from the final expression to obtain a 
somewhat simpler formula for the family of two-soliton potentials: 

Vi(x) = 12[4 + 9a 2 - a\ - 6 sech 2 x + 12 v /o^tanhx 
— (a,2 — 1)((<22 + 4) cosh2x + ^^/~a~2 sinh 2x)} 

x [a 2 -4 + (a 2 + 2)cosh2a; + 3a2sinh2x] _2 . (4.10) 

V 
-1 

-2 

-3 
-4 
-5 
-6 
-7 

-3-2-1 1 2 X 
Figure 1: Two soliton Poschl- Teller potential (dashed line) and its SUSY partner at ai = —4.1. 

For the potential difference to be regular one has to choose the factorization constant a 2 to 
be less than the ground state level of Vq, a 2 < —4, which gives complex values for the parameter 
d. Despite this, the potential difference is real and the solution ()4.8j) is also real (and nodeless) 
if it is set real at any point, i.e. ^(O) = 1- As an example we have plotted one of the potentials 
Vi = Vo + AV with AV as in ()4.9|) in the Fig. 1 (solid line) along with Vo (dashed line). 

5. The Lame-Ince potential 

A much more interesting example of second order Darboux displacements is given by the 
following result. 

Proposition 2 The periodic two-gap potential 

V (x) = 6p(x + uj') (5.1) 

where u' and uj are the imaginary and real half-periods of the Weierstrass function p, with g 2 
and g% as invariants, allows for second order Darboux displacements. The eigenvalues a\ >2 of 
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the transformation functions producing the displacement d are given by 



a l)2 = --[p{d)T\j92-^P 2 {d)}. (5.2) 

Proof. The proof consists simply in calculating the left hand side of the equation (|3.12|) . We 
obtain first some useful relations between our variables which are direct implications of the 
well-known properties of the Weierstrass functions (see e.g. [5]). From the differential equation 
for p one gets D" = SD, so T = (S'/D) 2 + 2S. The addition formula for p is 

V(u) - p'{v)~ 



p(u + v) + p(u) + p(v) = ~ 



(5.3) 



p(u) - p(v) 

By setting u = x + uj' + d, v = —x — to' and noting that p is even, this translates into 
(S'/D) 2 = + 4:p(d) and T = |5 + 4p(d). By taking the square root of the former of these 
relations and differentiating, we get 



'srv _ i 

D \ D ) ~ 3 



(5.4) 



and T = 8S'/3, D - (S'/D)' = 2D/3. Inserting these relations into (j3~T7j) yields 

p = 3(S'/D). (5.5) 

Now we are able to calculate the left hand side of (J3.12J) which gives 

Oi + 0,2 = —3p(d) . (5.6) 

This proves the first part of the statement. Now we shall find the eigenvalues of the 
transformation functions. For this purpose we shall use the equation (|3.1fij1 . Since [p'(x)] 2 = 
4p 3 (x) — g 2 p(x), we have the algebraic identity 

D 2 - 6^ + 3S 2 = 36g 2 , (5.7) 

and since p = 3S'/D, inserting it and ()5.6|) in (j3.16|) and taking into account (j5.7J) we obtain 

(a 1 -a 2 ) 2 = 9[g 2 -3p 2 (d)}. (5.8) 

This equality together with (J5.6|) proves the statement. One can check that formulas (J5.6j) 
and (|5.8|) . for the sum and difference of factorization constants, reduce to those for the specific 
case of the two-soliton potential ()4.5|) and ()4.6|) . respectively, when we take the half periods 
uo' = in jl and uj = oo. □ 

Some straightforward consequences of Proposition 2 are: 
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Corollary 2 Three of the five band edges for the potential , 

E x = 3e 3 , E v = 3e 2 , E 2 = 3d (5.9) 

correspond to d = uj, d = uj + uj' and d = uj' , where 

e 1 = p(uj), e 2 = p(uj + u/), e 3 = p(u/). (5.10) 

TTie lowest E = — \/3g 2 and the highest E 2 > = \/Sg 2 band edges are stationary points for the 
factorization constants ai and a 2 as functions of the displacement d. 

Corollary 3 The pair of (real) factorization constants which give rise to a real displacement 
by means of a second order Darboux transformation belong to the first finite forbidden band 
[Ei, Ei']. The real displacements so produced are in the interval (<i m i n , <i max ), where 

dmin = P~ X (vW 3 ) , d max = uj. (5.11) 

p~ l (t) is the function inverse to t = p(z) for which an appropriate sheet has to be chosen; 
d max = uj is realized for a 2 = Ei and a% = Ey , while rf min is realized for ai = a 2 = —a/3 g 2 /2. 

It is interesting to note that — are the positions of the band edges for the one-gap potential 
Vo = —2p(x + uj') and that in this case the range of the (first order) Darboux displacements is 
(0,oo). 

Proof. The first part of Corollary 2 follows immediately from Proposition 2 and known 
properties of the quantities given in Ref. [H] (vol. 2), when one uses d = uj, d = uj + uj' 
and d = uj'. The second one is a direct consequence of ()5.2j) and conditions a' 12 (d) = 0. 
Corollary 3 follows from the restriction on d to be such that g 2 — 3p 2 (d) be positive, since 
according to Proposition 2 it is inside of the square root. □ 

It is not difficult to see that for d = uj, the factorization constants are a 2 = E\ and a\ = Ey. 
This is the only possibility for the transformed potential to be displaced by half a real period 
of the Weierstrass p function. The possibility for the potential (|5.1|) to be displaced by a 
half-period was previously noticed in Here we indicate a possibility of displacing the 

argument of the potential (|5.1|) by a complex value and, in particular, by half an imaginary 
period. We remark that the values for the band edges are given in [2] for the Jacobi form of 
the Schrodinger equation, which we relate to the Weierstrass form in Appendix A. 

We are now able to integrate equations (|3.8|) - (|3.11|) and find analytic expressions for the 
Bloch functions, as it is established in the following proposition. 
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Proposition 3 The functions 



MxHW ^'V ) f^4) «<'-*«"'■ (5.12) 

4(5?) = ^\ (5.13) 
u 2 (x) 

are the Bloch eigenf unctions for the Hamiltonian with the eigenvalue 02, and 



a(x + uj') V 3C(d)x 



^>-""» U^^ J e (5 - 14) 

<M = (5.15) 
u x \x) 

are the Bloch eigenf unctions with the eigenvalue a±. Here a and ( are standard Weierstrass 
functions, 

xi = P~ l \ -\pQ + \yJ92 ~ 3po) - w ' 

> < (5.16) 

X2 = p~ l ( -\po -\y/92- 3poJ - <*>' 
po = p(d), p~ x is the function inverse to p and 

b = \((u , -x 2 )-±((u , -x 1 ). (5.17) 

The proof can be found in Appendix B. 

Remark 4 The solutions given in (\5.1 6 A ) and $5.13() are valid for any eigenvalue a 2 ■ By using 
(\5.0A ) one can find a\ as a function of a 2 : 



-a 2 ±\/9g 2 - 3al 

ai = y - • (5.18) 

Then with the aid of (j.5.61 ) d is expressed in terms of a 2 also. Hence, the Bloch functions 
<\5.1S\ ) and LH ) are determined only by the eigenvalue a 2 which can take any value and, thus, 
these functions generate the 1-parameter family of 2-SUSY invariant potentials mentioned in 
Corollary 1. With appropriate modifications the same is true for the functions < \5.14 ) and 

To illustrate this Remark we prove: 
Proposition 4 The potentials 

V 1 (x) = V (x)+AV(x), V (x) = 6p(x + u'), (5.19) 
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where 

[p(x + u' + d) + p(x + to')] 2 + 2p{x + cu' + d)p(x + uj') - g 2 /2 
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AV(x) 



p{x + uj' + d) + p{x + uj') + p(d) 
+ 3p(x + uj' + d) - 3p(x + uj') - p(x - xi) + p(x - x 2 ) 



(5.20) 



d 



ai+a.2 ' 



, X\, x 2 are defined as in (|<5.16|) . and a% in i\5. lt^i . form a one-parameter 
family of real-valued two-gap potentials isospectral with Vq. The role of the parameter is played 
by ci2 < Eq, where Eq = —y / 3g2 is the lowest band edge for Vq. 

Proof. The proof consists simply in calculating the right hand side of Eq. (jH.8j) and taking 
its derivative which will give us the potential difference (see (|2.4|) ). Using the expression (jB.2|) 
below for p in terms of the Weierstrass ( function and the addition formula (|5.3J) we first get 

(P - P) 2 



;(logp)' 



(5.21) 



2 V ° J/ P' + P' 

where p = p(x+uj') and p = p(x+uj' +d). To derive the final expression ()5.20|) we have used the 
formula (jB.8|) and the following properties of the Weierstrass functions: (p') 2 = 4p 3 — g 2 p — <?3, 
p" = 6p 2 — gijl and ('(x) = —p(x). The fact that these potentials have two-gaps follows from 
the fact that Darboux transformations preserve band structure (see e.g. ^1]). □ 

It is interesting to observe that different representatives of this family look like displaced 
copies of the initial potential (see Fig. 1 A) though Vq is not 1-SUSY invariant which is clearly 
seen from Fig. 1 B, where the left hand side of equation (J2.13|) is plotted. 



V 




-4.8 




-4.9 




-5 


\ : \ 




\: 1 


-5.1 


' A / 




• V/ 




-3 







A 



x 




B 



x 



Figure 2: (A) Lame-Ince potential (solid line) at uj = 1 and uj' — 2i and its SUSY partners. Dotted 
line corresponds to = —5 and dashed line to = —6; (B) illustrates that the left hand side of 
the equation 12.131 is not constant for d = 1.3. 
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V 

-5 
-5.2 
-5.4 
-5.6 

-6 -4 -2 2 4 X 

Figure 3: Potential corresponding to the transformation function u = + 0.5m \ at UJ = 1 and 
ui' = 2i and ct2 = —5.2 (solid line). Dashed line shows the initial Lame-Ince potential. 

Figure 3 shows a non-periodic potential obtained when a linear combination of the functions 
(J5.12j) and (J5.13j) is chosen as the transformation function. We stress that it possesses an 
energy level at —5.2 and the solution of the corresponding Schrodinger equation can easily be 
obtained by applying the first order Darboux transformation operator ()2.5|) to the functions 
(|5.12|) - (|5.15j) . We would also like to mention that the transformed potential tends to a displaced 
version of the initial potential at large values of |x|. This effect, first noticed in [3], illustrates 
a nonlocal deformation produced by a Darboux transformation. 

Finally, we would like to point out that the function (|5.12j) acquires a constant factor when 
its argument increases by the period 2u. This can be seen from properties of the Weierstrass 
a functions [5] that lead to 

u 2 (x + 2uS) = u 2 ( x ) ex P[C( a; )(3(i + x 2 — X\) — u(3((d) — 2b)} . 

From here it follows, in particular, that for a 2 = Se^, a\ = 3e 2 , one has d — u, X\ — u and 
x 2 = —2uj' which results in the condition u 2 (x + 2u) = u 2 (x) exp[2((u')uj — 2u'({^)} = —u 2 (x), 
where we have used the Legendre relation [S]. This result agrees with the fact that this function 
corresponds to the second band edge. 

In conclusion, we have characterized and investigated the interesting class of periodic 
potentials which are merely translated by a second order Darboux transformation. We have 
presented several examples, concentrating on the fascinating family of two-gap Lame-Ince 
potentials. The nonlinear differential equation which must be satisfied by all such potentials, 
which may be related to the class of Painleve equations, deserves further investigation. 
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Appendix A 

Here, to fix our notation, we relate the Weierstrass form of the Schrodinger equation to its 
Lame form (see also (2 , vol. 3). According to [5] one has 

2/ M ei-e 3 / , /\ , ( e 3 - e 2 )(e 3 - ei) 

sn (z,k) = —- , p{x + uj)=e 3 + (A.l) 

p(x) - e 3 p[x) - e 3 



or 



p(x + (J) = e 3 + (e 2 - e 3 )sn 2 (z, k) , (A. 2) 

where 

A; 2 = m= £^i ) z = (e 1 -e 3 ) 1 ' 2 x. (A.3) 
ei - e 3 

Now, after replacing the variable x by z in the Schrodinger equation with the potential (|5.1|) . 
one gets 

d 2 



dz 2 



+ (6k 2 sn 2 (z, k)-E) 



^ = 0, (A.4) 



where 



E = ^^ = — — + 2m + 2. (A.5) 
ei - e 3 ei - e 3 

Using this relation and Corollary 2 one recovers the known band edges for the Lame equation 
(see e.g. 0), 

~ E a — 6e 3 

E = — = 2m + 2 - 25, 

ei - e 3 

E 1 = m + 1, E v = Am + 1, (A.6) 
E 2 = m + 4, E 2 > = 2m + 2 + 25 . 
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In this Appendix we prove Proposition 3. In our notation the addition theorem for the 
Weierstrass ^-function (see [3]) 

1 p'(x + d) + p'(x) 



C(x)-((x + d) + ((d) (b.i; 



becomes 



Hence 



2 p(x + d) - p(x) 
^ = ((x + u' + d)-((x + u')-((d). 



P = 3^ = 6[C(x + oj' + d)- CO + J) - ((d)] . (B.2) 



Recalling that 

p 2 = 9 ( :_ ) =QS + 36p{d) (B.3) 



5 ,N 2 



D 



we find 



S = 6[p(z + a/) + p(z + w' + d)] . (B.4) 
Let us find now the expression for 1/p. 

The function S has second order poles at the same points as p{x + u') and p{x + u' + d), 
that is at x = —u' and x = —uj' — d. This means that S is a second order elliptic function 
and has two second order zeros. Here and below, for simplicity, we suppose that d falls in the 
interval (<i min , d max ) indicated in Corollary 3. 

Let us denote 

f(x) := p(x + u') + p(x + u>' + d) + p(d) . (B.5) 

Then p 2 = 36 / > for real values of p. Note that for d = u>, /(0) = f(u) = e\ + e<i + = 0. 
Hence, x = and x = w are the points of local minima for f(x). Therefore /'(0) = f'{oS) = 0, 
i.e. the zeros are of second order. Next, it is easy to see that f'(ou/2) = f'(3u/2) = and the 
points x = uj/2 and x = 3u/2 are the points of local maxima for real valued f(x). We find now 
the positions of these zeros. 

Let Xq be a minimum (or a zero-point) for f(x), that is f(x ) = and f f (xo) = 0. Let us 
abbreviate the notation by putting 

p := p(x + (J) , p := p(x + u' + d) , := . 
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Then p + p + p = and p' + p' = (with the evident notation p' = p'(x ), . . . ). Or else 
p' 2 — p' 2 = 0. Using the differential equation for p, p' 2 = Ap 3 — g 2 p — #3, one finds from here 
that 

(p - p)(Ap 2 + Ap 2 + App -g 2 ) = 0. 
Since p 7^ p this implies that 

= Ap 2 + Ap 2 + App -g 2 = Ap 2 + Ap 2 + App - g 2 

and 

p(x + u') = ~po ± -sjg~ 2 - 3p 2 Q . 



(B.6) 



Hence, the positions x\ and x 2 of the zeros of f(x) are given by (j5.16|) . In using these formulas 
it is necessary to keep in mind that p is a multisheet function. Note, that for d = uj one gets 

to + lu'. Note also that 



x 2 = and x 1 = u since y/g 2 - 3pg = e 2 - e 3 , p l (e 3 ) = uj', p 1 (e 2 ) 
since the equation p[z) = A defines z up to a sign and up to the periods, either the sum or the 
difference of X\ and x 2 modulo periods is equal to — d. 

The function p (x) = ±6 a/ f(x) has simple zeros x\ and x 2 and simple poles l/p(x). It is 
not difficult to find the residues of 1/p (x) at these points 

1 „ 1 1 



Res 



— Res 

X=X1 p [X ) x=x 2 P (X) Q^/ g2 _ 3p 2 

This implies that the residues of (a\ — a 2 )/p at the points X\ y2 are ±1/2. 

Now, using the known j3] decomposition of an elliptic function in terms of ( functions one 
finds 



(B.7) 



ai - a 2 
p 



^({x - ari) - ^({x - x 2 ) + b . 



(B.8) 



The constant b here may be calculated by the same formula at x = uj' where 1/p = 0. This 
gives ()5.17|) . 



Using the fact that 



ai - a 2 



EM. 

a(z) 



V 

and from (IB. 21) one finds 



log 



one obtains 

X\ 



a(x 



aix - x 2 



+ b 



p = 6 



log 



a(x + uj' + d) 
a(x + uj') 



6C(d) 



(B.9) 



(B.10) 
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With the help of the last two relations one gets from (J3.8)) 

[logu 2 = -[logpj — - — 
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(B.11) 



1 r, -1/ 3 

-[logp] +- 



log 



a(x + a/ + d) 



-few 4 



log 



cr(x — Xi) 
cr(x — X2) 



+ 6. 



cr(x + a;') 

It follows from here that (|5.12j) is a quadrature of the above formula. Similarly, the formulae 
()5.13|) - ([5.15p follow from f)5. 12j) and (j3.9|) — f)3.11j) . and the proposition is proved. 
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